A detailed analysis of structure growth in f{R) theories of gravity 
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We investigate the connection between dark energy and fourth order gravity by analyzing the 
behavior of scalar perturbations around a Friedmann-Robertson- Walker background. The evolution 
equations for scalar perturbation are derived using the covariant and gauge invariant approach and 
applied to two widely studied f{R) gravity models. The structure of the general fourth order per- 
turbation equations and the analysis of scalar perturbations lead to the discovery of a characteristic 
signature of fourth order gravity in the matter power spectrum, the details of which have not seen 
QQ ■ before in other works in this area. This could provide a crucial test for fourth order gravity on 

' cosmological scales. 
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■ In spite of all the efforts made so far, the problem of the nature of Dark Energy (DE) is still far from a completely 
satisfactory resolution. Among the many theoretical frameworks proposed, the idea of a geometrical origin of Dark 
Energy has recently received a great deal of attention. The main reason for this popularity can be found in the fact 
that these type of theories of gravity, which arc suggested by the low energylimit of very fundamental schemes P, [2] , 
lead to cosmologies which admit naturally a Dark Energy era [1, 0, [^, [a, 01 (and possibly even an inflationary one 
O ' [1, without the introduction of any additional cosmological fields. 

^ Much of the investigation performed up to now on the idea of Geometric Dark Energy has been focused on the 

so-called fourth order theories of gravity. In these theories the Hilbert-Einstein action is modified with terms that 
are at most of order four in the metric tensor. The features of fourth order gravity have been analyzed with different 
methods [ill . [T3 , [Tsl . [T^ . [Tsj and it has been shown that their cosmologies can give rise to a phase of accelerated 
expansion which is considered the footprint of Dark Energy. 

Although these results are very encouraging there are still some important open problems to be addressed. One of 
them is the analysis of the evolution of the linear perturbations and their comparison with observations. Over the 
\^ past year this problem has been studied by a number of authors, by (1) considering different ways of parameterizing 
' the non-Einstcin modifications of gravity or (2) by simplifying the underlying fourth-order perturbation equations 
I using a quasi-static approximation or a combination of (1) and (2) (igI. [TtI. llslTigt . 

' In a number of recent papers (20l . l2l| we derived the evolution equations for scalar and tensor perturbations of a 
00 subclass of fourth order theories of gravity characterized by an action which is a general analytic function of the Ricci 
. scalar. In our work we study the dynamics of linear scalar perturbations using the covariant and gauge invariant 
^ ■ approach developed for General Relativity (GR) in [H, [H, [H, [H, [l^, HI]. This approach has the advantage 
• i-H , of using perturbation variables with a clear geometrical and physical interpretation. Furthermore, we use a specific 
' recasting of the field equations that will make the development of the cosmological perturbation theory even more 
^ , transparent, allowing one to integrate the perturbation equations exactly for a given f{R) model without making any 
■ - - ■ additional approximations. 

The preliminary results obtained in [20| showed some interesting features. First of all the evolution of scalar 
perturbations is determined by a fourth order differential equation rather than a second order one. This implies that 
the evolution of the density fluctuations contains, in general, four modes rather that two and can give rise to a more 
complex evolution than the one of General Relativity (GR). Secondly, the perturbations arc found to depend on the 
scale for any equation of state for standard matter (while in GR the evolution of the dust perturbations are not 
scale dependent). This means that, for example, in this framework the evolution of super-horizon and sub-horizon 
perturbations is different. Third, and more surprisingly, we found that growth of large density fluctuations can occur 
also in backgrounds in which the expansion rate is increasing in time. This is in striking contrast with what one finds 
in GR and what one would naively expect, but at the same time suggests new ways to tackle the DE problem. 

The features mentioned above imply that the evolution of perturbations in this framework can be completely 
different from the one we are familiar with. Yet this does not necessarily mean that they are incompatible with 
observations. Rather, they are a sign of the fact that in dealing with these models one has to resist the temptation 
of using assumptions which work well in GR. In this paper, following in this spirit, we will analyze further what 
was found in 20] with the aim of achieving a clearer understanding of the physics of the matter dominated era in 
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fourth order gravity. In order to do this, we will rewrite the perturbation equations in a more physically meaningful 
way and will develop a series of tools which will make the analysis of the evolution of density perturbations easier to 
understand and to compare with GR. 

The paper is organized as follows. In section II we will give some basic equations and we will present briefly the 
covariant gauge invariant formalism we use to develop the perturbation theory. In section III, we give the background 
and the perturbation equations. In section IV we rewrite these equation in an interesting form allowing us to discuss 
their general structure. In section V we propose some useful tools to understand the behavior of the perturbations 
and compare it with what one obtains in General Relativity. In section VI we apply these tools to some simple specific 
examples. Section VII is dedicated to the conclusions. 

Unless otherwise specified, natural units {Ti = c = ks — SttG — 1) will be used throughout this paper, Latin indices 
run from to 3. The symbol V represents the usual covariant derivative and d corresponds to partial differentiation. 
We use the —,+,+,+ signature and the Ricmann tensor is defined by 



(1) 



where the W^bd are the Christoffel symbols (i.e. symmetric in the lower indices), defined by 



Wbd = -j^a"^ igbe,d + 9ed,b - 9bd,e) ■ 



(2) 



The Ricci tensor is obtained by contracting the first and the third indices 



Rab — g'^'^Rac 



(3) 



Finally the Hilbert-Einstein action in the presence of matter is given by 



A 



\R + Lr^ 



(4) 



II. GENERAL EQUATIONS FOR FOURTH ORDER GRAVITY. 



In four dimensional homogeneous and isotropic spacetimes i.e. Friedmann Lemaitre Robertson Walker (FLRW) 
universes, the most general action for fourth order gravity can be written as an analytic function of the Ricci scalar 
only: 



A= j d^X^[f{R)+Cm] , 



(5) 



where Cm represents the matter contribution. Varying the action with respect to the metric gives the generalization 
of the Einstein equations: 



1 



1 



fGab = f{Rab-:, 9abR = + ^9ab (/ ^ Rf) + ^b^af - Qab^ c^^ f , 



(6) 



where / = /(i?), /' = 



andT- = 



2 5{\F-9^^m) 



represents the stress energy tensor of standard matter. 



dr ' " ^9ab 
These equations reduce to the standard Einstein field equations when f{R) = R. It is crucial for our purposes to be 
able to write (El) in the form 



rpm I rpli rjitOt 

ab — -Lab ^ ^ab — ^ab i 



where Tj! 



r 



and 



T — — 

ab- J, 



^9ab if ~ Rf) + ^b^af - 9ab^c^'f 



(7) 



(8) 



represent two effective "fluids": the curvature "fluid" AS\ (associated with T^) and the effective matter "fluid" (as- 
sociated with T^) [3, [3,[23|- This step is important because it allows us to treat fourth order gravity as standard 
Einstein gravity plus two "effective" fluids. The details of the conservation properties of these effective fluids have 
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been given in [20'|. In particular, it has been shown that, no matter how the effective fluids behave, standard matter 
still follows the usual conservation equations = 0. 

The form ([7]) of the field equations allows us to use directly the covariant gauge invariant approach [l^, [H, [U, [H, 
m, [13, [13 in the same way presented in j20|. As usual the first step is to choose suitable frame, i.e., a 4- velocity Ua 
of an observer in spacetime. Following [20|, we will choose this frame to be the one comoving with standard matter, 
which is also called matter energy frame and will be indicated in the following by . We will also assume that in 
standard matter is a barotropic perfect fiuid with equation of state p = wp. Since the real observers are attached to 
galaxies and the galaxies follow the standard matter geodesies, this frame choice appears to be best motivated from 
a physical point of view. 

Once the frame has been chosen the derivation of the kinematical quantities can be obtained in a standard way [12] ■ 
In particular the derivative along the matter fiuid fiow lines is defined hy X — it^VX and the projected covariant 
derivative operator orthogonal to u° is given by V a = h^a^b- With these definitions we can define the key kinematic 
quantities of the cosmological model: the expansion Q, the shear Uab, the vorticity Uab and the acceleration aa — Ua- 
The general propagation equations for these kinematic variables in any spacetime correspond to the so called 1+3 
covariant equations [l^ which are given in Appendix |^ 

The definition of a frame Ua also allows us to obtain an irreducible decomposition of the stress energy momentum 
tensor. In a general frame and for a general tensor Tab one obtains: 



Tab = fJ-UaUb + phab + 2q(^aUb) + T^ab , 



(9) 



where /x andp are the energy density and isotropic pressure, qa is the energy flux (qa — q{a)) and TTab is the anisotropic 
pressure (iTab = 7r(ab>)- 

In this way, relative to w™, T*^* can be decomposed as 



T, 



^bc "a" 



la- 



rptot i^ab 
j^ab 

rjitot UC id 

^cd "■<a"-6> 



.R 

ab 7 



(10) 

(11) 



with 



/' 



/' 



a 



'ab 



"ab 
/' 



(12) 



Since we assume that standard matter is a perfect fluid in m™, g™ and are zero, so that the last two quantities 
above also vanish. 

The effective thermodynamical quantities for the curvature "fiuid" are 



'ab 



1 



-f"\/''RVaR-T:f R 
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/' 



f"Vi^a^b)R + rV(aR^b)R- labR 



lry'R+ 



(13) 

(14) 
(15) 
(16) 



The twice contracted Bianchi Identities lead to evolution equations for /i™, fi^ , q^ and are given in Appendix [XI 



III. LINEARIZED SCALAR PERTURBATIONS EQUATIONS 

Using the quantities defined above, and the equations given in Appendix |21 we are able to write both the evolution 
equations for the background and ones for scalar perturbations. As in [20| we will consider a background that is 
homogeneous and isotropic, i.e., a FLRW model. In this background the cosmological equations for a generic f{R) 
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read: 



3R 
2 



A™ + e (/i™ + p") - , 



0, 



(17) 

(18) 
(19) 



where R — 6K/S^ is the 3-Ricci scalar, K ~ 0, ±1 and 5 is the scale factor. The structure of these equations shows 
clearly that the effect of the introduction of higher order gravity on the background has a twofold nature. On one 
side, higher order gravity behaves like a an additional fluid in the model. On the other, it influences the way in which 
standard matter interacts gravitationally. 

Following [20j . we characterize scalar perturbations using the variables 



(20) 



the first three variables, which are borrowed from GR, represent the scalar fluctuations in the matter energy density, 
in the expansion rate (which is associated with /i) and in the spatial curvature. The last two represent the fluctuation 
of the Ricci scalar and its momentum. 

It is a relatively easy task to derive the propagation equations for these variables in a FLRW background. Their 
form is Appendix [Bl If we focus on the evolution of scalar part of these variables, which is associated with the 
spherically symmetric collapse terms, these equations become 



z = 



{l+w)Z , 



n 



Rf" 26 \ 

_ 1 i^r _ r Ai 
2 /' /' /' /' 

w ■ 

-R Am , 



(w-l)(3w + 2) fi 



2(ii;+ 1) 



RO 



/' 
RQ- 



2we2 

1 


+ 3w{n^ 






6{w + 1) 


f(3)' 
1 


n- 


w ~ 



/' 


w + 1 



A, 



/' 



3fJ 



f" - 
J V72 



/' 



1 



e + 2R 



(3) 



?ft-RZ 



1 f f(4) . . f(3) 



R 



{'3w — 1) /i 

3 r 
/^'^ R ^ 

7^ 

/" ■ 



3(w + 1) 



R 



A„ 



(21) 

(22) 
(23) 

(24) 



C 



18/"7^ 18A„ 



K 



C 



2{w-l)Q G/i^ 



1 



e/' 



6/" 



emff^^^ - f" (3/ - 2 (92 - 3/i^) /' + GRQf" 



e(/') 



l\2 



■n 



3(ii; + l) 



A, 



/' 

252/" 



/' 



252 [ef" - 3i?/(3) 



n 



together with the constraint 




Z - 2-^A,, 



2i?9^-^(/-2A. + 2i?9/" 



^5^-^V27^ = 0. (25) 



Note that this system is made up of four first order differential equation, which means that the evolution of every 
single perturbation variable is determined by a fourth order differential equation. This has a profound influence in 
the dynamics of the perturbations and makes them potentially very different from what one obtains in standard GR. 

Traditionally the analysis of the perturbation equations is simplified by using a harmonic decomposition. In the 
1+3 formalism this can be done by developing the scalar quantities defined above using the eigenfunctions of the 
Laplace-Beltrami operator (23 |: 



52 



(26) 
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where k = 2ttS/\ is the wavenumber and Q = 0. Developing ([20]) in terms of Q, (I2m25p reduce to 



Rf" 26 

T 



/' 



{w - l)(3u) + 2) /i 2wQ'^ + 3w(/z« + 3p^)"i 



2{w + l) f 



2 52 2/'/' /'/'^ 



6(w + 1) 



. f(3) 



sftC^-) = - ( 6 + 2i? 



w + 1 " 



(3w — 1) /i 



rfc2 
52 



1 f f(4) . . f(3) .. f{3) U 



/" 3(w + 1) 



R 



6i?e/7(3) _ 6A;V7' + /" (3/ - 2 (92 - 3Ai«) /' + 6i?e/' 



u; + 1 



e/" 



0(//)2 

4..2e_^^,^^2^^^,^_2^(er-3i?/(3; 



7^ 



'52 



3(w + 1) 



/' 



3/' 



(27) 

(28) 
(29) 

(30) 



/' 



3?+ 



(31) 



52 



/' 



y^2 



(32) 



which is a system of ordinary differential equations. This system takes a more manageable form if we reduce it to a 
pair of second order equations: 



1 



72 \„ i?r' 


AW _ 







fc2 
^52 



«;^-«;(3p^ + m'')- 



2wR&f" (3w2 _ 1) ^' 



/' 



/' 



A(fe) = 



(w + 1) 



2^f'-i + (f-'2t^ + 2Ref 



772 - 2^07r 



.2 

52 



3 6/' 6/' 



(3w- l)Ai+ 



1 + w 



/(3)7?2 + (p« + ^«)/' + l7?e/' 



+j?r^lAW-(^-w^^AW 



(33) 



(34) 



In the /(i?) = R case these equations reduce to the standard equations for the evolution of the scalar perturbations 
in GR: 



AW -{w-D eAW - 

7^W^(3^-l)MAW. 



fc2 



A^*") = 



(35) 



If one compares the system p3ll34p with the equations for the evolution of scalar perturbations for two interacting 
fluids in GR one notices that they have the same structure, i.e., one finds friction terms and source terms due to the 
interaction and the gravitation of the two effective fluids. It is then natural to ask ourselves if this analogy can be 
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useful to better understand the physics of these models. The answer is affirmative, but with some very important 
caveats. First of all a more correct way to draw this analogy would be to write the system of equations for 

and Afl — ^ and analyze their structure rather than using the ones above. Also, as stated in [23|, one has 

to be careful in remembering that we are dealing with effective fluids and, as such, they might violate some basic 
constraints that standard fluids usually follow (such as the energy conditions) or present subtleties in the definition 
of their comoving frame. 

However, in spite of these differences one can still use the coefficients of the (A™, A^) equations to obtain information 
about the interaction between standard matter and the curvature fluid. Unfortunately the length of this system makes 
it impossible for us to present it here (instead we will give their structure in the specific examples of Section VII). 
However, the coefficients of these equations are found to behave as a ratio of polynomials in the wavenumber and have 
a non trivial behavior in t. This kind of behavior is very different to what is found in a GR-two fluid system. In a 
photon-baryon system, for example, the dissipation terms grow as and behave like 1/t in time. This implies that, 
unlike Thompson scattering of the baryon-photon system, the effect of the interaction between matter and non-linear 
gravitation can influence large and small scales alike, depending on the structure of the action. Therefore, from the 
distribution of the structures in the observed sky one can deduce constraints on the nature of the theory of gravity. 
In fact, we will find that there is a specific spectral signature of fourth order gravity which is associated with these 
features. Another important difference with GR is that (|33II34I) are scale dependent for any value of the barotropic 
factor. This means that whatever the equation of state of standard matter, the perturbation solutions will always 
depend on the scale for which they are calculated, even in the special case of dust which in GR is associated to a scale 
invariant spectrum. 



IV. PERTURBATIONS AND COSMOLOGICAL PARAMETERS 



An interesting way of understanding the properties of the perturbations equations given in the previous section is 
to write them in terms of the cosmological parameters: 

9 = -9e-2|, ^7^ = -9^e-^ j^21^Q-\ (36) 

where q is the deceleration parameter, U,k the spatial curvature density parameter and j and s are the higher order 
kincmatical parameters jerk (or jolt) j and snap s (ssj . These quantities were devised in order to characterize the 
kinematics of a cosmological model in a way that is independent on any assumption on the dynamics and as we will 
see, they will be very useful for our purposes. 

Using g, J, s, and ^Ik defined in (|36p the Ricci scalar and its derivative can be rewritten as 

a2n _ / 2q 2\ 2 



R = Riq, nK,e') = (^y - 6f^A- + 3 j e^ (37) 

i^^ (1 + 1+417.-063, (38) 

V 27 27 27 3 9j ^ ^ 

In terms of these quantities the coefficients of the ((33l - [34|) above can be written in the form QP^F{q,j, s, ri^c, ki, Q^) 
where pi is a suitable integer associated with the dimension of the coefficient and ki — a~^Q^P^ are the wavenumber 
of the physical scales of the theory associated with the dimensional constants of dimension p2 in the action. In 
this way ([551 - become 

a(^> ^ e Ai^'^ +Be^ Ai^'^ c 7^('=) +ve ii'^^^ , (40) 
+ £ e Ttc^) 7^('=) = -g q^a^J^^ -ne^ A^^') , (41) 
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with 



2{j-q + 2rtK-2) f'Q^ 
9/' 



(42) 



B = 



5262 



2{w-l){j-q + 2nK-2)f"Q^ , 1 , , + 
9? + 3''(^ ~ + 



(43) 



{w + 1)/" fc2 (u; + 1) (2g2 + + 5 _ + 2)w + ji-q + '3w + 2) + (-4g + 6m; - 2)^^ + 2) /" 



9/(«; + l)2 



V 



(M; + l)(g(3K7 + 2)-j) 

292 {j -q + 2VIk -2) 4 (j - g + 2f7K - 2) 9*/' 
(^« + 1)9/" 
/' 



3 (j - g + 2nK - 2) /' 



(44) 
(45) 



ol(. _o 2(3g + j(g + 3) + g) \ 9(^^ + 1)/ (3M;g + g + 2»j, ~ 2) r 

3^'^ j-(Z + 2r!K-2 +7 2(j-g + 2r!K-2)94/" (j - g + 2r!^ - 2) 92/" ' ^ 

_^^_fc^_ 1 . 2_3^^_g_2-2f}-^ ■ 4(-j +9-20^ + 2)2/(4)94 ^ (3g + j(g + 3) + ,s)(2g - 3^« + 1) 



5*292 9^'' """" ^ " "-"Kj , , 9(-j + g - 217k + 2) 

(3(7 + jC? + 3) + s)2 3/(m; + 1) (3j + 9(9 + 5) + s - 2^1]^) , [2j' - (9(9^« + 7) + 2)j + 4(g + 2)^)2] /' 



9 (-j + 9-217^ + 2)2 4(-j + 9-2f7A' + 2)2/"92 6 (-j + g - 217^ + 2) 292/" 

[q {q^+q + s + 3(g(g + 5) + s)w - 18) - 2(.s + 4) + 2(-j + s + g(ll - Sqw) + 8)»k] f 

692 (_^- + ^ _ + 2) 2/" ' 

_ 2(w(3w - 4) + 1) (j - g + 2fiA- - 2) (3gw2 + 4f7A:w - 4w + g) /' (1 + 3w)/ 
^ 27(w + l) ^ 9(u> + l)92/" 692/" ' 



(47) 
(48) 



2(i.-l)(j-g + 2»j^-2) 

^ 9(u- + l) ' ^ ^ 

where / and its derivatives with respect to R are considered functions of R (g, r^x, 92) and ki. 

There are some general remarks that we can make at this point. First of aU, the fact that the perturbation equations 
can be shown to depend on the higher order cosmological parameters is a symptom of the fact that fourth order gravity 
is much more sensitive to the features of the background than GR (whose equations depend only on 9). Also, since 
these parameters can (at least in principle) be measured, we have a natural way to constraint both the dynamics 
of the background and the formation of structure. Note that this would not be possible in a standard scalar tensor 
theory of gravity, due to the fact that the background equations in that case remain second order. This suggest that, 
although in the background a scalar field might be able, in some cases, to emulate the behavior of fourth order gravity, 
this become more difficult at first order in perturbation theory [2^ [30| . 

Another important point concerns the meaning of the long and short wavelength limit. As it is clear from (j35p . in 

GR these limits can be defined by comparing the values of g-iQ2 and the matter term in the coefficient. However, 

looking at psp it is clear that the situation here is more delicate. For example, the short wavelengths regime cannot 
,2 7.2 

be defined as simply -g§Q2 ^ 1, but ^292 has to be bigger than all the other quantities appearing in the coefficients 

.2 

B, C and JF. The same reasoning holds for the long wavelengths: -^Jqi has to be smaller than all the other quantities 
appearing in the coefficients B, C and J- [i^. This effectively suggest the presence of a least three different regimes 
in the evolution of the perturbations. Firstly, the "deep super-horizon" regime in which k is effectively zero, an 
intermediate one (or two depending on the value of the barotropic factor w) which is determined by the details of the 
background and a "deep subchorionic" regime in which one has effectively k 00. 

The situation seems to become more complicated when the fourth order gravity action posses dimensional constants 
(like in the case f{R) = R + aR^ that we will consider later). Since these constants are associated with the scales 
at which the different contributions to the action become dominant, one would expect the introduction of additional 
scales into the theory, i.e., further possible evolution regimes for scalar perturbations. It turns out, however, that these 
additional constants lead only to changes in the power spectrum at intermediate scales, preserving scale invariance 
on large and small scales. In section VII we will see, using some examples, why this happens, and how these features 
can be used as a signature of fourth order gravity. 
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V. PROPERTIES OF THE SCALAR PERTURBATIONS AND COMPARISON WITH GENERAL 

RELATIVITY 

In order to understand the details of the evolution of scalar perturbation in a specific f{R) model we need to analyze 
the behavior of the solutions of the (|40j|4ip of that model. This can be done by examining their time dependence, but 
also by defining some characteristic quantities which help extract physical content from these solutions. One of these 
quantities is the power spectrum, P{k) of Am and TZ appearing in (|41l) . i.e., the variance of the amplitudes of their 
Fourier transform at a given value of k. In the case of (on which we will focus our analysis) P{k) is defined by 
the relation [39| : 

(A„.(ki)A,„(k2)> = P(fci)<5(ki + k2) , (50) 

where k^ are two wavevectors characterizing two Fourier components of the solutions of (|40ll4ip and P(ki) — P{ki) 
because of isotropy in the distribution of the perturbations. This quantity tells us how the fluctuations of matter 
depend on the wavenumber at a specific time and carries information about the amplitude of the perturbations (but 
not on their spatial structure). In GR the power spectrum on large scales is constant, while on small scales it is 
suppressed in comparison with the large scales j40l |. However, in the case of pure dust the matter fluctuations are 
scale invariant, so the power spectrum can be considered constant. This quantity is a powerful tool for comparing the 
predictions of the system p3ll34p with observations and is able to reveal a great deal of information on the physics of 
scalar perturbations. In addition, the analysis of the time variation of P{k) provides information on the ways in which 
the perturbations evolve in time on different scales. The fact that different scales evolve differently is a key feature 
of Geometric Dark Energy and if observed would allow one to differentiate between these models and ones based on 
standard Dark Energy. For our purpose we will normalize P{k) such that it is unity on super-horizon scales. This can 
then be scaled with current observations of the power spectrum on large scales (see ^31j] for the latest constraints). 

Another interesting way of exploring the properties of the equations presented above is to compare the features 
of their solutions with the corresponding results given in GR. This can be done by defining some suitable quantities 
whose value is associated with specific properties of the perturbation evolution. The first quantity is the ratio of 
the density perturbations in the modified theory with that of the corresponding result in the Einstein-de Sitter 
model, i.e., g = Am/^Eds ■ Since the evolution of scalar perturbations in the Einstein-de Sitter cosmology is 
AecIS — c+S + C-S~^^^, at late times it is approximately given by 5 « Am/S. The importance of the quantity g 
is due to the fact that g'^^ = 1 at all times. This means that calculating the quantity X = g — g^^ = 17 — 1 is 
a natural way of determining the deviation of the behavior of scalar perturbations in fourth order gravity from the 
corresponding result in GR. 

The second quantity is given hy Y ~ ^t\- This quantity is also measure of the deviation from GR, but has 

aln(S^) 

the additional advantage of measuring the behavior of scalar perturbations as a function of a at late times. This is 
particularly useful when one is forced to performed numerical integration of the equations above when exact solutions 
are not available. 

In what follows, we will use the quantities presented above to analyze the properties of two simple classes of fourth 
order theories of gravity. 



VI. EXAMPLES 



A. The case /(i?) = X-R" (i?"-gravity) 

The case /(i?) = X-R" also called sometimes i?"-gravity is characterized by the action 

L = ^[xR'' + Cm] , (51) 

and constitutes the simplest possible example of fourth order gravity. Its homogeneous and isotropic cosmologies 
have been studied in detail using the dynamical system approach [ill . 143 , l46l and the evolution of the large scale 
cosmological perturbations of a FLRW background has been investigated in 20] using the covariant gauge invariant 
approach. The results show some profound differences between this theory and GR. For example, in the transient 
Friedmann background a — i2ii/3(i+i«)^ .^g^g fg^j^j grow on long wavelength for almost all the values of the 
parameter 71, even when these values corresponded to backgrounds which undergo accelerated expansion (see Figure 
[T]for a plot of the real part of the exponents of the modes of the A^ solution). In this section we will continue the 
investigation of the evolution of scalar perturbations for this model [H^l focusing more specifically on the small scales. 
Substituting the form of f{R) in the general equations (|33ll34p we obtain the system (|CmC2p given in Appendix [Cl 
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FIG. 1: Plot against n of the real part of the long wavelength modes for i?"-gravity in the dust case (blue, red green and yellow 
lines) together with the GR modes (red and purple line). This graph is consistent with Figures (2] |3] and |4l i.e., when one of 
the four perturbation modes grows faster than GR X is bigger than zero and when they are all smaller than GR X is negative. 
Of course the cumulative effect of the combination of more than one growing mode has to be taken in account in comparing 
the two graphs. As one can see this produces a slight shift in the value of n associated with the change of sign in X and the 
dominance of one of the i?" growing mode over the GR growing mode. 



Providing the details of the background, the values of the parameter n, the barotropic factor w, the spatial curvature 
index K and the wavenumber k one is able to numerically integrate this system to obtain the behaviour of the matter 
fluctuations. We can then use the quantities described in the previous section to extract physical information about 
the evolution of density perturbations in this model. 

Figure [5] shows the behavior of the quantity X as a function of the time parameter r = log]^Q(S') on large scales 
{k « 0) for dust and different values of n. If we start with values of n close to 1+ the perturbations seems to first 
evolve faster than the growing mode of the GR-Einstein-de Sitter case and then at a slower rate which continues to 
decrease. For higher values of n, this first phase is absent and the rate is always slower than in the GR-Einstein-de 
Sitter case. For n > 3/2 the rate of growth of the perturbation is initially below the GR-Einstein-de Sitter mode, 
however, at late times the mode starts to grow at a much faster rate. For < n < 1 the situation is radically different: 
the rate of growth is always faster than GR. As n goes from zero to 1^ the growth rate seems to decrease, but at 
n ~ 0.85 this behavior changes and hy n = 0.90 the perturbations grow much faster than for any other value of n. 
This indicates the possible existence of an instability of the theory in this region of parameter space. 

These features are well summarized by Figure 3(a) which represents X at the present time r = 1 as function of n. 
When compared with Figure [1] one can see that the two figures are consistent. 

In Figure [3] we see the behavior of F as a function of t = logj^Q(5) for large scale perturbation. It is clear that the 
growth rate of the perturbations approaches asymptotically a constant value, which corresponds to the dominance of 
the fastest growing mode of the solutions of the system (|CmC2[) . Again, we found consistency between Figure [T] and 
both this plot and Figure |3(b)[ 

This confirms that X and Y provide a very useful source of information about the time evolution of scalar pertur- 
bations even when we do not have an analytic solution as it often happens in fourth order gravity. 

Let us then use the quantity X to analyze the growth rate of the matter fluctuation at different scales (see Figure 
[5]). One can see we see that the rate of growth decreases with respect to GR as one tends towards smaller scales. In 
fact, from the analysis of the time evolution of the power spectrum, we will see that this is associated to a dissipation 
of the perturbations at that range of scales. 

These results are consistent with the conclusions already drawn in [2^, i.e. that the dynamics of perturbation 
growth can be radically different from what happens in GR. This, in turn, implies that the time needed for linear 
structures to form and, as a consequence, the duration of the matter era, could be very different in this framework. 
Such features could prove to be an interesting method of placing constraints on the theory of gravity by using data 
from existing and future large scale structure surveys. 
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FIG. 3: Plots of X and F as a function of n evaluated at r = log]^Q(S') — 1 for _R"-gravity in the long wavelength limit. The 
horizontal lines X,Y — represent the value of these quantities in GR. 



It is also interesting to notice that the behavior of the perturbations for rt « 1 is also radically different from the GR 
case, i.e., even small deviations from n — 1 are able to produce completely different dynamics. One can understand 
such behavior since even small deviations from linearity in the action result in a change of order (from two to four) 
of the perturbation equations, so it is natural to expect significant differences in their solutions. 

As already mentioned the equations (|33ll34p and consequently (jCHIC2[) are similar in structure to the equations one 
would obtain in the classic multi-fluid GR case and these similarities can be used to infer the interaction properties 
of standard matter and the curvature fluid once the proper variables have been chosen (i.e. and A^j ). The full 
expressions of these terms in the background given above is too long to be reported in full. However it is instructive 
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(b)Plot of Y{t) for n = 1.55 



FIG. 4: Plot of Y{t) as a function of r = logi(,(S) for J?"-gravity with n > 1 (left) and < n< 1 (right). The red line Y = 
corresponds to GR. Note the consistency between this picture and Figure [1] 



to examine their structure. In the equation, the dissipation terms (i.e. the A'^ and the A^ terms) are given by 



Ca' = 



and 



t 



A2{n, w) 



9— 4ti / 

n,w) il + A5{n, w)k'^t 

Ei{n, w) 



t 



1 + E2{n,w)k^t 3(i+™)(i + ^3(n,w)A;2r ) 



respectively, while the source terms are 



and 



€ah 



Bi{n,w) 



Ci{n, w) 



1 + B2in,w)Pt ^oT^{l + B3{n,w)kH W^) 
1 + B4(n,-u;)fc2t^"W^(l + 85(71, w)kH'^~^i^) 



1 + C2{n, w)Pr~T»T+^{l + Csjn, w)kH^~'^^^^) 



(52) 



(53) 



(54) 



(55) 



where Ai, Bi, Ci, Ei are functions of n and w only. The structure of the coefficients of A'^ is, as expected, similar. 
For large scales (fc — > 0) we have 



Ca' 



Ai{n, w) [w + A2(n, w)] 



t 



Ei(n,w) Bi{n,w) Ci{n,w) 



t 



t2 



(56) 



which corresponds to the equations that one would obtain in GR with two coUisional fluids. For small scales {k ^ 00) 
we have 



Ca' --^^i^, Ca,-0, CA.-^i^ 



By,(n,w) 



B5{n,w) 



Ci{n,w) 



C5{n,w) 



(57) 



which again corresponds to the equations that one would obtains in the case of GR with two fluids, with the difference 
that this time, since Ca^ = 0, they are non-coUisional f5l!|. 
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(a)Plot of X{t) for n = 1.4 and different values of k (b)Plot of X{t) for n = 1.55 and different values of k 




(c)Plot of X{t) for n = 0.9 and different values of A; 

FIG. 5: Plot of X{t) as a function of r = logj^Q(S') for 7?"-gravity and k 0. Note how the growth of the perturbations is 
suppressed when k grows i.e. when smaller scales are considered. 



The form of these coefficients is very different to the ones obtained, for example, in the case of a GR baryon-photon 
system. As we mentioned in this last case the dissipation terms grow as k^. This means that, as expected, the effect of 
the interaction between the two fluids becomes more important when one considers smaller scales. In the above case 
the situation is different because the interaction scale of standard matter and curvature fluid is a non-trivial function 
of k which is peaked at a certain value of k. This means that the effect of the interaction is maximized around this 
specific scale. 

Let us now focus on the impact of these features on the matter power spectrum (Figures 6(a)[ |6(b)[ [7|). As 
mentioned in Section IIVI the k-structure of Equations (|33ll34p suggest that in fourth order gravity there exist at least 
three different growth regimes of the perturbations. This is confirmed by our results. In particular, in the case of 
dust we have three regimes for any values of the remaining parameters: (i) on very large scales the spectrum goes like 
GR i.e. it is scale invariant; (ii) as k becomes bigger the scale invariance is broken and oscillations in the spectrum 
appear; (iii) for even larger k the spectrum becomes again scale invariant. However, on these scales the spectrum can 
contain either an excess or deficit of power depending on the value of n. In particular for n « 1"*" small scales have 
more power than large scales, but, as one moves towards larger values of n, the small scale modes are suppressed. For 
< n < 1, instead the drop in power seems to decrease as one moves from n = towards 1~ and we see a sudden 
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increase for n « 1~. It is worth noting the case n « 0.8 for which there is basically no difference in power between 
large and small scales and there are no significant oscillations in the spectrum. 

Further indication of the link between the k structure of (|33ll34p and the different regimes of the matter power 
spectrum can be seen if one analyzes the power spectrum in a radiation dominated era. In this case the perturbations 
equations contain an additional term which is not present in the dust case. This means that one would expect four 
different regimes, rather than three. In Figure [8] we have plotted the matter power spectrum of i?"-gravity in the case 
of radiation and n = 10 (this value of n is chosen only for convenience and it does not have any physical motivation), 
and as expected one can recognize four different regimes. Note that we obtain the same k scaling as in GR on small 
scales. 

Finally further information on the dynamics of the matter perturbations can be obtained examining the time 
evolution of the power spectrum. In Figure [H] we give the power spectrum for n = 1.4 at different times. One can see 
that, as the universe expands, the small scale part of the spectrum is more and more suppressed and oscillations start 
to form. This is in agreement with what one finds from the analysis of X and suggests that in this model small scale 
perturbations tend to be dissipated in time. On the other hand the large scales do not seem to be evolving, which 
might appear in contrast with what mentioned above. However this is a byproduct of the normalization: for clarity 
we have normalized the spectrum in such a way that every curve has the same power in long wavelength limit. 

The features of the spectrum that we have derived can be then interpreted in terms of the interaction between the 
curvature fluid and standard matter. On very large and very small scales, the coefhcients (j52ll55p become independent 
from k so that the evolution of the perturbations does not change with the scale and the power spectrum is scale 
invariant. On intermediate scales the interaction between the two fluids is maximized and the curvature fluid acts as 
a relativistic component whose pressure is responsible for the oscillations and the d issip ation of the small scale per- 
turbations in the same way in which the photons operate in a baryon-photon system |52l | . The result is a considerable 
loss of power for a relatively small variation of the parameter n. For example, in the case n = 1.4 the difference in 
power between he two scale invariant parts of the spectrum for n = 1.1 is of one order of magnitude while for n = 1.6 
is about ten orders of magnitude. 

Probably the most important consequence of the form of the spectrum presented above is the fact that the effect 
of these type of fourth order corrections is evident only for a special range of scales, while the rest of the spectrum 
has the same k dependence of GR (but different amplitude). This implies that we have a spectrum that both satisfies 
the requirement for scale invariance and has distinct features that one could in principle detect, by combining future 
Cosmic Microwave Background (CMB) and large scale surveys (LSS) [32l. i33|. 

B. The case f{R) ^ R + ai?" 

We will now consider a theory with the following action, 

L = [R + aK" + Cm] ■ (58) 

This theory h as g ained much popularity as a fourth order gravity model within the context of both inflation and dark 
energy BESP- 

Unlike JJ^-gravity, (j58p includes explicitly the Hilbert-Einstein term, so that one can consider it as the result of an 
additive correction to GR. This also means that this model introduces an explicit physical scale which is determined 
by the relative magnitude of the two terms in the action, making it the simplest fourth order gravity theory with an 
additional scale for the gravitational interaction. 

Some work on estimating the parameter a for the specific case of n = 2 was completed during the eighties. This 
work was mainly based on arguments pertaining to black hole physics and quantum gravity (see for example _44]). 
In the following we will take a, which in our units is the ratio between the coupling constant of the fourth order 
corrections, to be positive definite. Of course, we expect this model to behave much in the same way as the model 
discussed in the previous section for a — > cx), and to recover GR for a = 0. This also means that any new feature in 
this model will emerge for intermediate values of the coupling. 

The f{R) = R + aR^ model has been analyzed at the level of the background using many different approaches 
(see for example d [13, SHE!), but probably some of the most interesting results for cosmology have been found 
using the dynamical system approach [la . |4^ . l47l | . The dynamical systems analysis proved that this class of models 
has, like i?"-gravity, an unstable fixed point associated with the Friedmann-like solution a — i2n/3(i+iu) Substituting 
the form of the action ([55]) , the system of perturbation equations (|33ll34p takes the form (jDmD2p given in Appendix 
Id] and one can analyzed the evolution of the scalar perturbations in this background [ssj . 

In spite of the fact that both this theory and _R"-gravity have a Friedmann-like solution, dealing with the pertur- 
bation equations for this background in /(i?) — R + ai?" is considerably more complicated than in i?"-gravity. As 
a consequence, one is unable to find exact solutions for the perturbation equations, even when one adopts the long 
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(a)Plot of the power spectrum at r = 1 for /{"-gravity and ra > 1. Note that the spectrum is composed of three 
parts corresponding to three different evolution regimes for the perturbations. 
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(b)Plot of the Power spectrum as a function of k for i?"-gravity at t = 1 
and < n < 1. The spectra for n = 0.6 and n = 0.5 only approach the 

scale invariant plateau at extremely high k when compared to the other 
curves. Note the behavior of the spectrum for n fa 0.8, differently from all 
the other cases, there is basically no loss of power in the spectrum at large 



FIG. 6: Power spectra for _R"-gravity 
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FIG. 7: Detail of the plot of the power spectrum for 7?"-gravity at r = 1. As expected, for these values of n we find the 
presence of the three regimes mentioned in the text. Note also that for n ~ 1"'' the small scales are characterized by an excess 
of power. Such features is compatible with the results found in [l^. [itI. [l^ . 
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FIG. 8: The power spectrum for i?"-gravity at r = 1 in the case of radiation for GR and ii"-gravity with n — 10. As expected 
in this last plot we find four different regimes instead of the three of the dust case: a first regime for — > which is scale 
invariant; a second and a third regime which correspond to the two different slopes between k « 10~^ and k ~ 10^^^ and a 
fourth regime which has the same slope of the GR plot. 
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FIG. 9: Evolution the Power spectrum for _R"-gravity for n = 1.4. The spectrum has been normalized in such a way that the 
curves coincide at large scales. Note how, as time passes, small scale perturbations are dissipated and oscillation appear. 



wavelength limit. To progress one must numerically integrate the system of perturbation equations and extract useful 
information relating to the perturbation dynamics using the quantities X, Y and the power spectrum. 

If we look at Figure [TO] we notice that for a > 1 the dynamics of the large scale perturbations in this model are 
similar to those found in i?"-gravity but contain some small differences. For example, from Figure 10(a) one sees that 
when a — 10 and n « 1+ the growth rate has some oscillations at small r followed by the decay typical of i?"-gravity. 
However, at smaller values of a the situation changes. In particular, the curves that did not show oscillations at small 
T are characterized by the onset of late time oscillations. Because of that the value of n at which the growth rate 
starts to increase at late time is different. These features reveal the important dynamical differences between this 
model, i?"-gravity and GR. Such differences will be even more evident when we will examine the time evolution of 
the power spectrum. 

The analysis of the behavior of the quantity Y on large scales also reveal similarities between this model and i?"- 
gravity. In particular, one sees that, again, for large values of a (i.e. a = 10) there is very little difference between the 
previous example and this model, but when a becomes smaller the solution takes more time to saturate and oscillates 
with a wider amplitude. In addition, the asymptotic values of Y that correspond to the dominant perturbation mode 
changes with a and n. 

Let us now consider the evolution of the perturbations on smaller scales. The behavior of X at different scales 
(Figures [T^ reveals changes in the growth rate. However for different values of a and n these changes are not always 
associated with suppression of the growth rate like in the case of i?"-gravity. Figure |12(d)| gives an example of these 
differences. 

When one calculates the structure of the coefficients for the Am equation in the system (A„i, A/j) one obtains 



Ca' - 



and 



Ai{n,w,A) 
t 



A2{n,w,A) 



Ei{n,w,A) 



(59) 



t 



1 + E2{n,w,A)k^t^~^^^{l + E3{n,w,A)k^i^~^^^) 



(60) 



(a)Plot of X{t) for f(R) = R + aR" for large scales 
(fc = 0), n > 1, dust and a = 10. Note the similarity 
with the plots of Figure [2] 
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(c)Detail of the plot |10(a)| for 1 < ra < 1.05, The 

oscillations in the curves are the feature that 
differentiate f{R) = R + aR" from ij"-gravity. 
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(e)Plot of X(t) for f{R) = R + aR" for large scales 
(k = 0), n > 1, dust and a = 0.1. 



(b)Plot of X(t) for f{R) = R + aR" for large scales 
(fc = 0), n > 1, dust and a = 1. Differently from the 
plot in Figure [lO(a)| in this case the curve n = 1.55 
reaches a maximum and then decreases. 
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(d)Detail of the plot |10(b)| for 1 < n < 1.05, The 

oscillations in the curves are the feature that 
differentiate f{R) = R + aR" from fj"-gravity. 




-2 -1 
Lo.E (T) 



11=1 n=1.05 11=1.15 n=1.25 ii=1.45 n=1.55 



(f)Plot of X(t) for f(R) = R + aR" for large scales 
(k = 0), n > 1, dust and a = 0.01. 



FIG. 10: Plots of X{t) for large scales in ii + ai?"-gravity 
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(a)Plot of y (r) for /(_R) = R + aR" for large scales (fc = 0), n > 1, (b)Plot of Y{t) for f(R) = R + aR" for large scales (fc = 0), n > 1, 
dust and a = 10. dust and = 1. 
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(c)Plot of y (t) for f{R) = R + aR" for large scales (fc = 0), n > 1, (d)Plot of Y{t) for f{R) = R + aR" for large scales (A; = 0), n > 1, 
dust and a = 0.1 dust and a = 0.01. 



FIG. 11: Plots of Y{r) for large scales in R + ai?"-gravity 



respectively, while the source terms are 



Bi{n,w,A) 
i2 



1 + B4{n, w, A)fc2r"3n+^(l + 35(71, w, A)fc2r"3(T+^o )_ 



(61) 
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(c)Plot of X{t) for n = 1.4, a = 0.01 and various (d)Plot of X(t) for n = 1.55, a = 0.01 and various 

values of k values of k 

FIG. 12: Plot of X{t) as a function of r = logj^o(S') for R + ai?"-gravity and k 0. Note how the growth of the perturbations 
is suppressed when k grows and when a decreases. 



and 

1 + C2{n, w, A)kH^-^(^{l + Csin, w, A)kH^-^(^)] 

4 4 ) (OZj 

1 + C4(n,ii;, A)fc2t^"3<T+^(l + C5(n,u;, A)fc2t^"wT^)J 

where Ai, Bi, Ci, Ei are functions of n and w and A ~ a t^"^^. Although the time behavior of these coefficients 
is always different from GR (unless of course a = 0), their k behavior at a fixed time is similar to the one found in 
the case f{R) = i?". This means that this model also behaves like a two fluids system for large and small k and that 
the deviation from scale invariance occurs only in a specific k interval, which this time is determined by both the 
parameters a and n. 

These features characterize the spectra shown in Figures [131 As one can clearly see, these plots resemble the ones 
we have derived in the previous example. There are, of course, differences in the position of the oscillations and the 
amount of the power drop, but one finds again three different regimes in the case of dust and two of them (fc — > and 
A; ^ oo ) correspond to scale invariance. Particularly interesting is the fact that in principle the values of a and n 



Ci(n, w. A) 

i2 
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can be fine tuned in such a way to obtain a spectrum in which the small scales have the same power as the large ones. 
In a situation like this most of the spectrum would be scale invariant and all the deviations would be concentrated 
around a specific scale. 

The time evolution of these spectra also reveals some interesting insight into the dynamics of the matter fluctuations. 
As usual for large values of a the evolution is very similar to the one obtained for i?"-gravity as it is shown in Figure 



14(a) However, when the value of the coupling changes the behavior of the perturbations can change dramatically. 
An example is given in figure [l"4(b)| in which the evolution for the power spectrum of the model (n = 1.4, a = 00.1) 
in which the small scale perturbations are first dissipated and successively start to grow again. This means that in 
principle one could choose n and a, such that for example the small scale perturbation grow at different rates at 
different times. This property could be useful in the resolution of open problems in GR structure formation, like the 
cosmological dark matter or the excess of red galaxies. 

In conclusion, in spite of all the differences in the dynamics of perturbations, the power spectrum in this class of 
model seems to preserve most of the main structure of the one in i?"-gravity. This implies that all the considerations 
made in the previous section concerning the physical mechanisms behind the form of the spectrum can be made also 
in this case. This result was not expected and suggest that we might have encountered a characteristic signature of 
/(i?)-gravity which would be crucial to investigate the validity of these models. 



VII. DISCUSSION AND CONCLUSION 



In this paper we have discussed the properties of scalar perturbations in fourth order gravity and described some 
useful methods for extracting physical information from the complex equations that govern their evolution. These 
tools are devised in such a way to gain as clear an understanding as possible of the behavior of scalar perturbations 
on all scales and to facilitate a direct comparison of these results with the corresponding ones in GR. Two simple 
models: f{R) = -R" and f{R) = R + aR" were analyzed in detail because of their simplicity and because we have a 
relatively good understanding of their background via the dynamical systems approach. 

The results obtained show profound differences in the dynamics of the perturbations compared to what occurs in 
GR: in these models we find the growth rate of the perturbations is in general different from the GR one and always 
dependent on the scale. This implies, for example, that depending on the value of the parameters the perturbation 
can grow or dissipate (or both) at different rates, with obvious consequences for the global cosmic history. 

The fourth order system of differential equations that governs the behavior of the scalar perturbations of these 
models, in particular when written in terms of Am and A^, has a structure that resembles closely one of a two fluids 
GR model. Although the A,„ and system of equations is enormously more complicated than (|33II34|) . one can still 
use it to get an idea of the nature of the interaction between the non-Einstein part of the gravitational interaction 
and standard matter. 

In section VII we found that the wavenumber structure of these coefficients is such that for very large or very small 
k they become scale invariant. This implies, in turn, that the matter power spectrum is scale invariant for fc — > 0,oo 
and can present some characteristic features on scales that depend on the different parameters of the model. Another 
way of interpreting the form of the spectrum without necessarily using the curvature fluid idea is to interpret TZ and 

as being associated with the propagation of the scalar degree of freedom of the theory, or a scalar gravitational 
mode, whose interaction with matter is able to emulate the effect of a relativistic component, which, like photons in 
GR, induces power loss in the oscillations that appear in the spectrum. 

The picture that emerges from our results is that in our examples, fourth order gravity influences deeply the 
structure formation process, but the modiflcations are only detectable around a speciflc value of k. Everywhere else in 
k space the results are very close (although dynamically different) to GR. This is particularly interesting since it means 
that for a suitable choice of values for the parameters, the oscillations can be positioned beyond the observational 
boundary of currently available data. Hence our results seems to imply that not only could these models be compatible 
with the observed matter power spectrum, but that we also have a systematic way of constraining their parameters 
using data coming from a range of scales, for example by combining CMB and LSS data jsil. [ssj. 

It is also worth commenting briefly about the compatibility of our analysis with the existing literature. For example 
in ^1^ a class of models which is very similar to the one we analyzed in Section IVI Bl is considered. Using a speciflc 
background the authors showed that the matter power spectrum is characterized by an excess of power at small scales 
when the theory is very close to ACDM. A similar result is found in [l7| in which classes of models are considered 
which contain additive corrections to the Hilbert -Einstein action (i.e. they have the form R + g{R) ). These general 
corrections are parameterized by a quantity B, which measures the deviation from GR. Instead in (l8| . generic 
modifications are considered (both scale-dependent and scale- independent), using a more flexible parameterization. 
In the sub-case of scale dependent modifications (like in the f{R) case) examples are given in which one finds once 
again excess of power on small scales in the power spectrum. Also one should note that in the same sub-case examples 
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(a)Plot of P{k) for n > 1, a = 10 and various values of k 



(b)Plot of P{k) for n > 1, a = 1 and various values of k 



Log [P(k)l 




Log |P(k)l 




I n=l n=Ll 11=1.2 n=L3 il=L4 ii=L5 n=1.6 n=L7| | 0=1 11=L1 ii=L2 il=L3 ii=1.4 n=L5 11=L6 ii=L7 

(c)Plot of P(fc) for n > 1, cj = 0.1 and various values of k (d)Plot of P{k) for n > I, a = 0.01 and various values of k 

FIG. 13: Plot of the Power spectrum as a function of k for R + Q_R"-gravity at r = 1 for n > 1. 



were given when one finds a deficit in power. 

Interestingly, both the examples considered in this paper exhibit excess power on small scales for n — > 1^ in the case 
of iJ^-gravity and for n 1+ and a <C 1 in the case R + aR^\ This indicates that, for small corrections to GR, excess 
power at large scale seems to be a generic feature in /(i?)-gravity. However there are situations in which one might 
want to analyze theories which are not necessarily close to GR. For example in [s^, [13, [IHl a fit of i?"-gravity with the 
data coming from the rotation curves of galaxies and supernovae type la leads to values of n in the range [1.7, 3.5]. 
For these values of n, both our examples indicate a loss of power at small scales, which means that we are provided 
with an opportunity to rule out this model by testing it against available data. Although the presence of an excess or 
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(a)Plot of P{k) for n = 1.4, a = 10 and evaluated at 
various values of r 



(b)Plot of P{k) for n = 1.4, a = 0.01 and evaluated at 
various values of r 



FIG. 14: The time evolution of the Power spectrum in _R + a_R"-gravity for n — 1.4 and differing values of a. Note the 
drastically differing vertical scales in the plots. Additionally, note the rise in power at small scales at late times in the case of 
n = 1.4 and a = 0.01. 



deficit of power on small scales seems closely related to the value of the specific parameters of the model itself, there 
are some indications that other features of the spectrum found in our examples are indeed general. For example, the 
k structure of the general equations (|40II41|) suggests that the evolution of perturbations in a generic f{R) theory 
presents at least three different regimes. Also looking at the derivation of the (Am,Afl) equations (which can be 
performed in general, provided a sufficiently large amount of paper and time) one realizes that the k dependence of the 
dissipation and source coefficients we have found in our examples is expected to be common to any f{R) Lagrangian 
because it originates from the V^V^ terms in the perturbation equations. 

We end by commenting that these results tog ether with the dynamical systems analysis of the background cosmo- 
logical history presented in other papers [Til . |46| provides a unified and consistent approach to the combined study of 
FLRW observational constraints and a complete analysis of linear structure growth in the context of f{R) gravity. 
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APPENDIX A: GENERAL PROPAGATION AND CONSTRAINT EQUATIONS OF THE 1+3 

COVARIANT FORMALISM. 



Expansion propagation (generalized Raychaudhuri equation): 

e + + a„(,<7»'' - 2a;„a;« - V»«„ + + + 3p™) = -^{n"" + 3p^) . (Al) 
Vorticity propagation: 

+ I + ^curluo - cTflbW^ = . (A2) 

Shear propagation: 

^{ab) + fQcTab + Eab - ^ {aUb) + CTciaO'b}'' + ^{a^b) " ^{aH) = ^T^ab ■ (■'^3) 

Gravito-electric propagation: 

Et^^b) + QEab - CmlHab + +p"')<Jab - 2u%d(a-H"6)'' " StJeia^h)' + ^''e^d(aEb)'^ 

= -slM +P )(^ab - 2^{ab) - 2^ {alb) " 6^^ab " 2^ {a^b)c " 2^ ^c{a^b)d ■ (^^j 

Gravito-magnetic propagation: 

H(ab) + QHab + CUrlEab — 3(7c{aHb}'^ + W^£cd(a-f^6) + '^^'^^ cd{aEb)'^ 

= ^curlTT^, - lu^^q^^ + \cj\aeb)c(ld ■ (A5) 

Vorticity constraint: 

V"a;„ - ii^iJa = . (A6) 

Shear constraint: 

VV„6 - CUrlWa - fVaO + 2[w, U]a = -Q^ • (A7) 

Gravito-magnetic constraint: 

curl (Tab + V{a^b) - Hab + '^U{a'^b) = • (A8) 

Gravito-electric divergence: 

V^Kb - I V„/i" - [a, HI + SHabOj'' = \alq^ - | [w, qX - 5 V^-tt^, + i VaM"" - l^Qa ■ (A9) 
Gravito-magnetic divergence: 

V'-Fab - ir +Pn^a + [a,i?]a - 3^6^^ = -icurlgf + (^^ +p^)u;a - W.^X - ^7r>^ . (AlO) 
Standard Matter Conservation 

/i" = - e (/i" + p") , (All) 

V V" = -(/x" + P^) ii" . (A12) 



Curvature fluid Conservation 

leq : cons2M« + V>f = - 6 (/x« + p«) - 2 («°gf ) - (a<*7rf„) + ^ , (A13) 

Zeg : cons3gf„> + V„p« + VV^, = - | e 5^ - a^" - (/.^ +p'')ua-u' tt^, - ^ '^6 «f + m'" ^^^^ (AM) 
As usual angle brackets applied to a vector denote the projection of this vector on the tangent 3-spaces 

V^a) = ha'Vb . (A15) 
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Instead when applied to a tensor they denote the projected, symmetric and trace free part of this object 

W^{a6> = [hia'h/ - ^h^'^hab] W^d ■ (A16) 

The spatial curl of a vector and a tensor is 

{X)'' = e"''" WbXc (X)''^ = e"''^" ^cX% (A17) 

respectively, where Cabc = u'^rjabcd is the spatial volume. Finally uja — \£a"^^bc and the covariant tensor commutator 
is 

[W,Z]a=eacdW\Z''\ 

The 1+3 equations above are completely equivalent to the Einstein equation and govern the dynamics of the matter 
and gravitational fields in fourth order gravity. As we will see the new source terms in their R.H.S. will modify the 
evolution of the perturbations in a non-trivial way. The standard GR equations are obtained by setting f{R) = R 
which corresponds to setting all these sources to zero. 



APPENDIX B: THE COVARIANT GAUGE INVARIANT GRADIENT EQUATIONS 



In the following we give, for completeness, the equations for the evolution of the gradient variables 
S 



Za = S^aQ, Ca^SVaR, Ua = SV aR , ^ SV aR ■ 



They read 



( Rf" 28 \ 



3(w - 1) (3ui + 2) ^ 2wQ^ + 3w{fi^ + 3p^) 2w k 



6{w + l) /' 



2 ^ /' 52 2 /' /' ^ /' /' ^ V /' 



6{w + l) 

/(3) 



152 



&f" 



2/' 



RQ 



f" ~9 



W -\- I f 



w + l 



Rv:\ 



(Bl) 
(B2) 

(B3) 
(B4) 



5Ra = 



. f(3) 

- ( e + 2i?^ ) SR, - RZa - 



(3w — 1) /i 



/' 



/" w + l j" 6(W + 1) 



3 , \ ^ 



3 /" /' 



— (w + 1) RIS — 

6 /' 2 ^ V 3 /' 



f -~ 9 
7^a + V27^a , 



/(3) 

/" 



(B5) 



together with the constraint 



Ca 2Rf" 



■ • f(3) f" / . k 

2RQjr - 772 - 2m + 2i?e/" - 4-^ 



7^a + ^^3f^a-^v^7^, = 0. 



/' 



The propagation equation for the variable C is 



(B6) 



Ca 



_ 6fc2 / / 



5*2 e 



Ca 



Qj^ _ 2(3c^)9 \ _ er 

e 3(c. + i)y' e/' 



6i?7"2 (2 (92 - 3m^) /' - 3/) /" 6i?7^^^ 



/'2 

4cj52e 



e/'- 



/' 



3(c. + l) " / 



2^2^/, 252 (^er - 37?/(3) 



3/' 



7^a 



Tea 



(B7) 



These equations were already given in [20|, but contain some typos that we have been corrected in the equations 
above. 
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APPENDIX C: PERTURBATION EQUATIONS IN THE CASE f{R) = xR" 
The general perturbation equations in the case f{R) = X-R" read 

All'J +A&Al^'>-Be^ Al^'> = cn^''^ + V ii^''^ , (ci) 
^t*^) +S& Tt'*^) - 92 7^('=) = -g e^A^^) -hq^ A w , (C2) 



where 



3(9-90^ + 1) 3' ^^-^^ 

wfc^ 2 (3(n - 2)w^ + (n - l)w - n + 2) (g - 90^ + 1) 

^ ^202 + 3;^^ 

(n - 1) (3^2 + w-l){j + q + ISQk - 2) ((n - 3)n + 2) (Sw^ +w-l){j + q + ISnK - 2) ^ 



q-9nK + l 3(g- gf^if + 1)2 

(n - 1) (3w;2 +w-l){{q-8)q + s + 18{q - 3)nK + 6) 



3(9-9f2^ + l) 



3(n- l)(w + l) 
2 ((J -917^ + 1) 



(2(n - 3)n + 3) (g - 9f2if + 1) (n(3n - 7) + 4) (j + q+ IS^k - 2) 



£'292 3(n-l)n 3(n - 1) (g - 9n/f + 1) 



(n - 2)(n -l){j + q + ISQk - 2) 2 (n - 1) ((g - 8)q + s + 18{q - 3)17x + 6) 



3 (g - 9f2K + 1) 2 3 (g - 90^ + 1) 



(C4) 



(C5) 



3(n-l)(w + l) 

^=2(,-90k + 1)' ^""'^ 

2{n-2){j+q+18QK-2) 

3(«-90x + l) ' ^^^> 

fc2 (n-2)(2n + 3(n-l)«;)(g-90if + 1) (-3w + n(3w + 2) - 4) (j + g + 180^ - 2) 

~ 5292 9(n- l)n 6(9-90^ + 1) 

(2 - n)(n(3M; + 2) - 3(w + 2)) {]+q + 180^ - 2) 2 {3w + 4- n{3w + 2)) ((g - 8)q + s + 18{q - 3)0/^ + 6,)^„, 

TSTZ — — r~r\ r^^) 



18(5-917^ + 1)2 18(g-9fi^ + l) 

2(n - 2)(3«; -l){q- Q^k + 1) ^ 3(«;(3w + A)-l){j + q+ 180a: - 2) 



{n — l)n w + 1 

(n - 2)(w(3w + 4) - 1) (j + g + 180^^ - 2) 2 _ (w(3w + 4) - 1) {{q -8)q + s + 18{q - 3)0^ + 6) 
(w + 1) (g - gOif + 1) w + l 



(C9) 



(CIO) 

9(«; + l) ^ 

APPENDIX D: PERTURBATION EQUATIONS IN THE CASE f{R) = R + aR" 

The general perturbation equations in the case f{R) = R + aR^ read 

A(^) + y^l 9 A(^) - S 92 A(^) = C7^('^) + 1? 9-1 Tte^' , (Dl) 
^^^'^^ + f 9 Tt^^) - 92 7^('=) = -e? 9''aW - H 9^ A(^) , (D2) 
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where 



4(j 



K 



2) 



3 (3fc„ + 4 (g - 9f7K + 1)) 



(D3) 



B = 



2 {{3 - 9w^) + 2w {q - QQk + I)) {q - 9flK + 1) A {3w^ + w - l) {j + q + ISflx - 2) 
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9kaW 
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2w (j + q + 180k - 2) 
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— 3w 
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2g2 4r>Ag 16g ^ 2s ,„ ^ 4 

— — 4ii&r H 
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